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ABSTRACT 

From radial velocities determined in high signal-to-noise digital spectra, we report the 
discovery that the brightest component of the binary system HD 150136 is of spectral 
type 03. We also present the first double- lined orbital solution for this binary. Our 
radial velocities confirm the previously published spectroscopic orbital period of 2.6 
days. Hen absorptions appear double at quadratures, but single lines of Nv and Niv 
visible in our spectra define a radial velocity orbit of higher semi-amplitude for the 
primary component than do the He n lines. From our orbital analysis, we obtain min- 
imum masses for the binary components of 27 and ISM©. The neutral He absorptions 
apparently do not follow the orbital motion of any of the binary components, thus 
they most probably arise in a third star in the system. 

Key words: stars: binaries, spectroscopic — stars: individual (HD 150136) — stars: 
O-type— 



1 INTRODUCTION 

HD 150136 (a(2000) = 16h41m24s; <5(2000)=-48°45'47" ; 
V=5.5), together with its close companion HD 150135, are 
the brightest stars of the open cluster NGC 6193 in the nu- 
cleus of the Ara OBI association. Their spectra have been 
classified by Walborn (1972) as 05III:n(f) and 06.5V((f)), 
respectively. HD 150136 appears as a double lined binary in 
the compilation of binary stars of spectral type O (Garmany 
et al. 1980), with a period of 2.7 days and mass ratio of 1.8, 
but the orbit has not been published yet. Arnal et al. (1988) 
in their search for short period binaries in the open cluster 
NGC 6193, published a single-lined orbit with a period of 
2.6 days. 

HD 150136 is also the brightest star of a visual multi- 
ple system with 6 components within a radius of about 20 
arcseconds. The closest component to HD 150136 is a 9th 
magnitude star (component B of the system) at 1.6 arcsec- 
onds (cf. Mason et al. 1998). HD 150135 (V=6.7, component 
C) is the second brightest star in this group at a separation 
of 9.6 arcseconds from HD 150136. 

In this paper, we report the first double-lined orbital 
solution for HD 150136 based on medium resolution digi- 
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tal high S/N spectral images. Furthermore, our spectra and 
radial velocities suggest that the primary star has a very 
early spectral type 03. At the distance of AraOBl, about 
1300pc (cf. Herbst & Havlen 1977), this star then becomes 
the nearest star of this early spectral type. 



2 OBSERVATIONS 

We have obtained 19 digital spectral images of HD 150136 
with the REOSC Cassegrain spectrograph attached to the 
2.1m reflector at the Complejo Astronomico El Leoncito 
(CASLEO 1 ) in San Juan, Argentina. The spectra were ob- 
served in May and June 2004. A TEK 1024 x 1024 pixels 
CCD was used as detector. The reciprocal dispersion of the 
spectra is 1.8A / pixel. Exposure times were between 1.5 
and 5 minutes, resulting in spectra of signal-to-noise ratio 
between 100 and 150. The wavelength range covered in our 
spectral images is from 3800 to 5500A . Comparison lamp 
spectra of CuAr were observed immediately after or before 
at the same telescope position as the stellar images. Flat- 
field and bias frames were also observed each night. Onedi- 

1 CASLEO is operated under agreement between CONICET, SE- 
CYT, and the National Universities of La Plata, Cordoba and San 
Juan, Argentina. 
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mensional spectra were extracted from the two-dimensional 
spectral images using IRAF routines. 



Table 2. Preliminary Circular Orbital Parameters for 
HD 150136. 



3 RESULTS AND THEIR DISCUSSION 
3.1 The spectrum of HD 150136 

Along with the characteristics of the previously determined 
spectral type of 05III:n(f), i.e. strong absorption lines of 
H and Hen, fainter absorptions of He I, and Niii emission 
at 4634-40A our spectra of HD 150136 show also absorp- 
tion lines of Nv 4603-19 A and Niv 5203 A as well as the 
emission line of N IV 4058 A . These lines are the distinctive 
ones observed in earlier 02-3 type spectra, (cf. Walborn et 
al. 2002; Gamen & Niemela 2002). One of our spectra of 
HD 150136 is illustrated in figure 1. The simultaneous pres- 
ence of He I and Nv absorptions in the spectrum clearly 
pointed to a composite spectrum and to the possibility that 
the brighter binary component is of earlier spectral type 
than 05. We therefore decided to study the radial velocities 
of the lines in the spectrum of HD 150136. 



3.2 The radial velocity orbit 

Radial velocities were determined fitting gaussian profiles 
to the spectral lines within the IRAF routine SPLOT. He II 
absorptions, and to a lesser extent also those of Hydrogen 
Balmer lines, have the appearance of double lines of unequal 
components in several of our spectra. We did not use the H 
absorption lines in our radial velocity study. Radial velocities 
for He II lines were determined using the deblending routine 
of SPLOT. The absorption lines of N V and N IV, as well as 
the Niv emission, appear single in all of our spectra, and 
we used the radial velocities of these lines to determine an 
approximate orbital period. The journal of our radial veloc- 
ity observations is presented in Table 1, and spectra corre- 
sponding to approximately opposite orbital phases (0=0.21 
and 0=0.71) are depicted in Figure 2, where the lines of 
both components can be appreciated. 

We introduced the values of Nv and Niv radial ve- 
locities in the period search routine of Lafler & Kinmann 
(1965). As can be seen from the values listed in Table 1, 
the radial velocities from our spectra show large variations 
from one night to another, but smaller variations between 
observations obtained during the same night, thus indicat- 
ing a binary period of a few days. We therefore searched for 
periods between 1 and 5 days. The best period we found was 
2.65 days, almost identical to the values found previously, 
namely 2.7 days and 2.6 days (cf. Garmany et al. 1980, Arnal 
et al. 1988). 

The value of the period we found was then introduced 
as an initial value to the routine for defining the orbital 
elements of the binary. To this end we used an improved 
version 2 of the program originally published by Bertiau & 
Grobben (1969). The radial velocity orbit of the primary 
star was derived using the mean of the radial velocities of 
Nv 4603-19A and Niv 5203 A absorptions and the Niv 
4058A emission line, since these lines likely originate only in 

2 available upon request from 

|ftp://li len.fcaglp.unlp.edu.ar/pub/fede/gbart-Q. 1-41. tar, gz 
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the brighter component. The orbit of the secondary compo- 
nent was determined from the radial velocities of the Hell 
absorptions observed when the components appeared most 
separated. To a first approximation, the radial velocity orbit 
appears to be circular, as we obtained an orbital eccentricity 
of 0.03 ± 0.02. Therefore we have fitted circular orbits to our 
radial velocity observations. The orbital elements are listed 
in Table 2. These orbital elements should be considered as 
preliminary values, because due to the rather limited wave- 
length resolution of our spectra, the separation of the spec- 
tral lines of the secondary component may not be adequate. 
The upper panel in Figure 2 illustrates the double-lined ra- 
dial velocity orbit of HD 150136. 

The minimum masses that we find for the binary 
components are moderately high (cf. Table 2). However, 
HD 150136 is not known as an eclipsing binary, and from 
Hipparcos photometry Marchenko et al. (1998) find a dis- 
persion of 0.03 mag, with no trace of periodicity. Thus we 
may expect the orbital inclination to be rather low. An in- 
clination of 50 ° would result in masses of about 60Mq and 
4OM0 for the 03 and 06 components, respectively. 

The radial velocities of He I absorption lines apparently 
do not follow the orbital motion of any of the binary com- 
ponents. However, in a few of our spectra, these absorptions 
show a fainter component to the red when the secondary 
of the binary system has its maximum positive velocity. We 
therefore assume that the He I absorptions mainly originate 
in a third star in the system, and thus HD 150136 is a multi- 
ple star similar to Sk-67degl8 in the Large Magellanic Cloud 
(Niemela, Seggewiss & Moffat, 2000). This fact is also sup- 
ported by the lower amplitude of the radial velocity vari- 
ations of the Hell absorptions of the primary component 
(138km s _1 ), as compared with the radial velocity varia- 
tions of N v and Niv lines (217km s _1 ). 

The middle panel and the lower panel in Figure 3 illus- 
trate the radial velocity variations of the stronger compo- 
nents of Hell and of Hel absorptions, respectively, phased 
with the binary period. 

A lower amplitude of Hen absorption radial velocity 
variations probably arises because despite of the possible 
line blending effects, these lines are also blended with those 
of the third star in the system. As noted by Andersen (1975) 
for the diffuse He I lines, pair-blending effects are appreciable 
on lines which show extensive wings, even if their cores ap- 
pear well defined and resolved. This effect probably also af- 
fects the He 11 absorptions. We have tried to correct for pair- 
blending of Hell absorptions in our spectra of HD 150136 
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Figure 1. Continuum rectified spectrum of HD 150136 obtained at CASLEO in 2004, May. Absorption lines identified below the 
continuum arc H.5A4102, H7A434O, H/3A4862, Hen AA 4200, 4541, 4686, 5411, He 1 AA4471, 4922, NvAA4604, 4620, and NivA5403. The 
emission lines identified above the continuum are N IV A4058 and N III AA4634-4640-4642. 



using the multiple profile deblending routines within IRAF, 
which fit gaussian, lorentzian, or voigt line profiles with a 
linear background. Obviously, to verify the success of this 
procedure in separating the lines of the secondary compo- 
nent, as well as to ascertain the contribution of the third star 
to the stronger components of the He 11 absorptions, an im- 
proved orbital solution based on high resolution observations 
is needed. However, the high masses of the binary compo- 
nents for an assumed orbital inclination of 50° together with 
the lack of photometric variations (see above), suggest that 
our preliminary values of the semiamplitudes of the radial 
velocity variations in the HD 150136 binary system may not 
be very far from the true values, unless the masses of the 
binary components are much higher than those determined 
for other stars of similar spectral types (e.g. Massey et al. 
2002). 



3.3 The spectral types of the binary components 

As noted in section 3.1, the presence of high ionization Nv 
and Niv lines in the spectrum of the primary component 
points to an early 02-3 spectral type. The emission and ab- 
sorption complex of N ill and He 11 at ~ 4634-4690 A changes 
its shape with orbital phase. Emission of Hell 4686 A is ap- 
parent when the primary is in front of the system. The spec- 
tral classification criteria for 02-4 stars put forward by Wal- 
born et al. (2002) compare the N IV emission at 4058 A with 
the N in emission at 4634-40 A . In our spectra of HD 150136 
this criterion would indicate a spectral type 03.5. The lu- 
minosity class is difficult to infer, but assuming that the 
Hell 4686 A PCyg type profile observed when the primary 
is in front of the system arises in this component, the pri- 
mary would appear to be an 03.5 If* type star. 

We have estimated the spectral type of the secondary 
component in the HD 150136 binary system from a com- 
parison of the spectral lines of He I and Hell 4471/4541 A 



4 V.S. Niemela & R.C. Gamen 

Table 1. Journal of spectral observations of HD 150136 
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Phases were calculated according to T = 2 453 171.2+2.662 E 



in our spectra when the secondary component has its maxi- 
mum positive radial velocity. This is illustrated in Figure 4, 
where these lines corresponding to the secondary are de- 
noted with dotted vertical lines. In this Figure we note that 
the He II 454lA line corresponding to the secondary is some- 
what stronger than He I 4471 A , indicating a spectral type 
earlier than 07, where the relation of these lines becomes 
unity. Together with the fact that Hen 4686A line, which is 
a luminosity indicator in O type stars, appears in absorption 
(cf. Fig. 2), the spectral type of the secondary component 
appears to be 06 V. 



4 SUMMARY AND CONCLUSIONS 

From an analysis of recent high signal-to-noise spectra of 
HD 150136 we find the following: 

(i) Radial velocity variations of spectral lines with a pe- 
riod of 2.662 days are present in our spectra of HD 150136, 
confirming the previously found values of the binary period 
(Garmany et al. 1980, Arnal et al. 1988). 

(ii) Radial velocity variations of the highest ionization 
lines, namely N v and N IV, define an orbital motion of rather 
higher amplitude for the primary star than the absorptions 
of Hen, 217±4 km s _1 vs. 138±5 km s" 1 , respectively. 

(iii) The radial velocities of the main components of He I 
absorptions do not follow the orbital motion of any of the bi- 
nary components. Therefore HD 150136 appears to harbour 
at least another luminous OB type component, in addition 
to the 2.6 day binary system. 

(iv) The primary component of this binary system is of 
spectral type 03, noticeably ealier than the spectral classifi- 
cation found in the literature, namely 05III:n(f), which ob- 



viously referred to the composite of all three spectra present, 
at a single phase. We estimate a spectral type for the sec- 
ondary component as 06 V. 

(v) Previously reported small light variations do not show 
periodicity, thus the orbital inclination is probably low. An 
inclination of 50 ° would produce masses of 60 and 40 Mq 
for the 03 and 06 components, respectively. 

HD 150136 is one more example of the increasing num- 
ber of O stars discovered to belong to tight multiple star 
systems, such as e.g. CPD -59° 2636 and CPD -59° 2603 in 
the open cluster Trumpler 16 in the Carina Nebula. The for- 
mer system consists of an 07 V+08 V short period (3.62d) 
binary and a fainter 09 V component (Albacete Colombo 
et al. 2002); and the latter system similarly of a short- 
period (2.15d) eclipsing 07 V+09.5 V binary bound to a 
B0. 2 IV star (Rauw et al. 2001). Another relevant example 
is HD 167971, a triple system with a short period (3.32d) 
eclipsing binary and a more distant star of spectral type 
08 lb, which is the most luminous component of this triple 
system (Leitherer et al. 1987). 

Close triple stars are also found among the Wolf-Rayet 
stars, thought to be evolved descendants of O type stars. 
Notorious examples of these are 9 Mus (= HD 113904) 
and HD 5980. 9 Mus consists of a WC+(0?) short period 
(18d) binary with a wider component which is a more lumi- 
nous late O type supergiant ( cf. Moffat & Seggewiss 1977). 
HD 5980 resides in the brightest Hn region of the Small 
Magellanic Cloud, and contains a Luminous Blue Variable 
(LBV) which erupted in 1994 ( Barba et al. 1995). The sys- 
tem consists of two emission line stars in an eclipsing binary 
of short period (19d) and a line of sight O type companion 
(cf. Niemela et al. 1999). 




Figure 2. Spectra of HD 150136 during opposite orbital phases (<p=0.2l and </>=0.71) are depicted showing the spectral lines of the two 
components. Lines of the primary are denoted with vertical solid lines and those of the secondary component with dotted lines. 



Mason et al. (1998) in their astrometric/spectroscopic 
survey of O stars brighter than V~ 8, find a large fraction 
of primary components in close visual systems to be short 
period spectroscopic binaries. These "hard" binaries in close 
triple and multiple systems are clues to our understanding 
of massive star formation processes, as well as the formation 
of runaway stars by gravitational encounters. 
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